Entanglement between an auto-ionization system and a neighbor atom 
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I. INTRODUCTION 

Ionization is a process in which an electron is trans- 
ferred from its bound discrete state into a continuum of 
free states, e.g., by interacting with an optical field. In 
a stationary optical field, an electron at an atom grad- 
ually leaves its bound state and moves into an ionized 
free state This process is irreversible. It can be uti- 
lized for the generation of entangled electron states that 
are stable in time. The timc-dcpcndcnt entanglement 
among bound electrons can easily be generated in re- 
versible interactions (Coulomb interaction, dipole-dipole 
interaction) . The irreversible ionization can subsequently 
'freeze' it and provide this way the stability in time. We 
demonstrate this approach on the simplest model of two 
atoms, one of which allows the electron ionization. 

It is well known that the process of ionization is 
strongly influenced by the presence of additional discrete 
excited states (auto-ionization levels). They consider- 
ably modify the long-time photoelectron ionization spec- 
tra (for an extended list of references, see, e.g. @4a|)- 
There even might occur Fano zeros 6-8] in the spectra 
of isolated auto-ionization systems due to the mutual in- 
terference of different ionization paths. The interaction 
of an auto-ionization system with neighbor atoms leads 
to the presence of dynamical zeros [9M13I] that occur peri- 
odically in time. Ionization spectra contain useful infor- 
mation about bound states of an atom and that is why 
they have widely been studied experimentally 14] . Auto- 
ionization systems have also been found useful as media 
exhibiting electromagnetically-induced transparency and 
slowing down the propagating light [l5j . The ionization 
process is also sensitive to quantum properties of the op- 
tical field • 

Here, we consider two atoms in a stationary optical 



field that moves electrons from their ground states into 
excited or ionized states. Electrons in their excited states 
mutually interact by the dipole-dipole interaction [l7| . 
This creates quantum correlations (entanglement) be- 
tween two electrons. Whereas one electron remains in 
a bound state, the second one is allowed to be ionized. 
We pay attention both to the temporal entanglement for- 
mation [l8|, [l9[ and its long-time limit. The quadratic 
negativity of a bipartite system generalized to the con- 
tinuum of states is used to quantify the entanglement. 
We show that highly entangled states can be reached in 
a wide area of parameters characterizing the system of 
two atoms. 

The paper is organized as follows. A semiclassical 
model of the system under consideration is described in 
Sec. II together with its the most general solution. The 
formula for negativity as a measure of entanglement in a 
bipartite system with the continuum of states is derived 
in Sec. Ill and compared with quantum discord. The 
spectral density of quadratic negativity is introduced in 
Sec. IV. The dynamics of entanglement as well as its long- 
time limit are discussed in Sec. V. The spectral entan- 
glement is analyzed in Sec. VI. Conclusions are drawn in 
Sec. VII. Appendix A is devoted to an alternative deriva- 
tion of the formula for negativity. 



II. SEMICLASSICAL MODEL OF OPTICAL 
EXCITATION OF AN AUTO-IONIZATION ATOM 
INTERACTING WITH A NEIGHBOR ATOM 

We consider an atom b with one auto-ionizing discrete 
level that interacts with a neighbor two-level atom a by 
the dipole-dipole interaction (for the scheme, see Fig.[IJ. 
Both atoms are excited by a stationary optical field. This 
composite system can be described by the Hamiltonian 
H, 
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FIG. 1: Sketch of an auto-ionization system b interacting with 
a two- level atom a. Ground states are denoted as |0} a and 
1 0) b whereas symbols |l) a , |l)f>, and \E) stand for the excited 
states. Dipole moments jj, a , Mti an d M describe the appro- 
priate interactions. The excited discrete state at atom a (b) 
has energy E a (Eb), whereas energies E characterize excited 
states \E) of the continuum. Symbol V describes the Coulomb 
configurational coupling between the excited states at atom b. 
Constants J a b and J emerge from the dipole-dipole interaction 
between the atoms a and b; oll is the pumping amplitude. 



Here, the Hamiltonian -ff a _; characterizes the auto- 
ionization atom: 



ff a _i = E h \\) hh (\\ + j dEE\E)(E\ 

+ JdE [V|£) b (l|+H.c.] 

+ \p b a L exp(-iExi)|l}6b(0|+H.c.] 

+ J dE [[ia L exp(-iE L t)\E) b (Q\+K.c.}. (2) 

Energy Eh means the energy difference between the 
ground state 1 0) b and the excited discrete state |1)& of 
atom b. Similarly, energy E stands for the energy dif- 
ference between the state \E) in the continuum and the 
ground state |0)&. The Coulomb configurational coupling 
between the excited states of atom b is described by V. 
The dipole moments between the ground state |0) & of 
atom b and its excited states are denoted as \i and /ih- 
The stationary optical field with its amplitude cll oscil- 
lates at frequency El- We assume % = 1. 

The Hamiltonian -fft-a of the neighbor two-level atom 
a introduced in Eq. ([1} takes on the form: 

fft-a = E a \l) aa (l\ + [»a exp(- iE L t)\ l) aa (0\ + H.C.] , 

(3) 

where E a means the energy difference between the 
ground state |0) a and the excited state |l) a ; /i a stands 
for the dipole moment. 

The Hamiltonian -fftrans in Eq. (TTJ) characterizes the 
dipole-dipole interaction between electrons at atoms a 



(Ja b |l>6b(0||0) aa (l|+H.C.) 

d£ [J|£) 6 <0| |0) ao (l| + H.c.] . (4) 



In Eq. (Q}, J a h (J) quantifies the dipole-dipole interaction 
that leads to the excitation from the ground state |0)b 
into the state 1 1) & {\E)) of atom b at the cost of the decay 
of atom a from the excited state | l) a into the ground state 
|0) o . 

Following the approach of Ref. [Ill ], a state vector 
\if>)(i) of the system at time t can be decomposed as 

IV'X*) = C d o(t)|0)a|0)6 + C 1 o(t)|l).|0) 6 

+ C Dl (t)|0)«|l) 6 + Cli(t)|l) o |l)6 

dEd (E,t)\0) a \E) 
dEdx{E,i)\l) a \E) (5) 



using time-dependent coefficients coo, cqi, cio, en, d${E), 
and di(E). 

These coefficients satisfy a system of differential equa- 
tions which can be conveniently written in the matrix 
form: 

d f 

—c(t) = -iAc(t)-i / dEBd(E,t), 

at J 



dt 



d(E,t) 



iB^c(t) -iK(E)d(E,t) 



(6) 



and 



c(t) 



coo(i) 
cio(i) 
coi(i) 
cn(t) 



d(E,t) 



d (E,t) 
di(E,t) 



(7) 



The matrices A, B, and K introduced in Eq. ^ are 
time-independent provided that a basis rotated at the 
pump-field frequency El is used: 











H a a L A a J* b 
H-bOLL Jab A & 

ii b a L li a a L A a + A b 



Vb a L 
V>L 



B 



K(E) = 



4 

/'- 


a i 


o - 


J* 


H*a* L 


V* 










v* 




' E 


-E L 



fJ-aC*L 



t^a^L 

E-E L + 



(8) 
(9) 

(10) 



Here A a = E a — E L and A b = E b — E L stand for the fre- 
quency detunings of discrete excited states with respect 
to the pump-field frequency. 
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Contrary to the solution of the model equations found 
in [ll| we adopt here the most general approach based on 
algebraic decomposition of dynamical matrices and solu- 
tion of the corresponding Sylvester equation. We first 
neglect threshold effects in the ionization, eliminate con- 
tinuum coefficients d(E) in Eq. (|6]), and introduce a new 
matrix M: 



M = A inBB^. 



(11) 



The matrix M describes the dynamics of only discrete 
states that is governed by the vector c. We denote eigen- 
values of the matrix H(E) as E ~ £1 and E — £ 2 and 
eigenvalues of the matrix M as Amj, j — f , 2, 3, 4. These 
eigenvalues occur in the matrix decompositions of matri- 
ces K and M: 



K(E) 
M 



4 

E 

3=1 



A-m-,-M,- 



(£-&)K 2 , (12) 
(13) 



The basis matrices Ki and K2 can be obtained from 
the following equations: 



K! + K 2 = I 2 , 

(£-&)Ki + (£-&)K 2 = K(E). 

The eigenvalues £1 and £2 are given as follows: 

_ A a ±5£ 
£1,2 - El 7. > 



(14) 



S£ = ^/Al+4\fi a a L \ 



(15) 



Symbol 5£ means the frequency of Rabi oscillations of 
the two-level atom a. 

Similarly, the basis matrices Mj, j — 1, 2, 3, 4, arise as 
the solution of the following equations: 

4 

E M , = 14, 



k = 1,2,3. 



(16) 



In Eqs. (O and ([IB]), I 2 and I 4 are 2 x 2 and 4x4 unit 
matrices, respectively. 

After the introduction of matrix M in Eq. (|lip. the 
solution of Eqs. (HJ for the vector c can be written in the 
very simple form: 



c(f) = exp(-iMt)c(0); 



(17) 



c(0) is the vector of initial conditions. 

On the other hand, a newly introduced matrix T{E) 
(of dimension 2x4) obtained as the solution to the 
Sylvester equation [20| 



K(E)T(E) - T(E)M = B* 



(18) 



is useful for expressing the solution of Eqs. (j6j) for the 
continuum of states described by the vector d(E). On 
using the matrix decompositions written in Eqs. (|12[) and 
(Ti"3"|) , the solution of the Sylvester equation (fl"8"|) can be 
expressed as follows: 



fc=i 2=1 



1 



A 



-KjB^M 



Mj 



J ■ 



(19) 



The components of amplitude spectrum of an ionized 
electron at atom b are given by the coefficients in the 
vector d(E,t). They can be written in the most general 
form 



d{E,t) = (exp[-iK(E)t]T(E) ~ T(E)exp[-iMt 



c(0) 
(20) 

depending on the initial conditions. We have assumed 
that d (E, 0) = d\ (E, 0) = in the derivation of Eq. ([20) . 

As the interaction processes between the discrete states 
and the continuum of states are irreversible, the eigen- 
values of matrix M are complex with negative imaginary 
parts. As a consequence, the expression in Eq. (|20p for 
the amplitude spectral components d simplifies in the 
long-time limit: 



d lt (E,t) = exp[-iK(E)t]T(E)c{0). 



(21) 



III. NEGATIVITY OF A BIPARTITE SYSTEM 
IN DISCRETE AND CONTINUOUS HILBERT 
SPACES 

We need to quantify the amount of entanglement be- 
tween the two-level atom a and the auto-ionization atom 
b that has a continuous spectrum. The philosophy 
based on declinations of the partially-transposed sta- 
tistical operators of entangled states from the positive- 
semidefinite partially-transposed statistical operators of 
separable states [2ll |22| has been found fruitful here and 
has resulted in the definition of negativity. 

Following the approach by Hill and Wooters [2l|, we 
write a matrix P of the statistical operator describing an 
electron at atom a and a (fully) ionized electron at atom 



b in a given time T [dj {E) 



dj{E,T),j 



0,1 



d (E)d* (E>) d (E)dl(E') 
di(£K(£') d x {E)dl{E') 



(22) 



We note that the frequencies E and E' in Eq. (f22|) are 
considered as continuous indices of the matrix P. 

A partially-transposed matrix P Ta transposed with re- 
spect to the indices of two-level atom a is obtained after 
the exchange of sub-matrices in the upper-left and lower- 
right corners of the matrix P in Eq. (|22jl: 



>Ta 



d (E)d* (E') d l {E)dl{E') 
d (E)dt(E') d 1 (E)d\{E l ) 



(23) 



In order to determine negativity N , we need to find 
the eigenvalues A of matrix P Ta first. An eigenvalue A 
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together with its eigenvector (uq(E), ui(E)) fulfil the fol- 
lowing system of equations with a continuous index E: 



d (E) J dE'd* (E')u (E') 

+ d 1 (E) J dE'd* (E')ui(E') = Xu [E), 

d (E) J dE'd{(E')u (E') 

+ d 1 (E) J dE' d\{E')ui{E') = X Ul (E). 

(24) 

Integrals in Eqs. (|24|) give the coefficients a,j k of the de- 
composition of eigenvector functions Uj (E) in the basis 
dj{E): 



a jk = / dE'd*(E')u k {E'), j, k = 0, 1. 



(25) 



Using the coefficients aj k defined in Eq. (j2"5"j) . the equa- 
tions in (|2"4"|) can be rewritten as follows: 

d (E)a 00 + di(E)a i = Xu (E), 

d {E)a 10 + d l {E)a ll = \ Ul (E). (26) 

The projection of equations in Eq. (j2"rJ)) onto the basis 
vectors dj (E) results in a system of four algebraic equa- 
tions for the coefficients ajk determining the eigenvector 
(u (E), Ul (E)): 



boo b i 

b Q0 b i 

bio hi 

b w hi 





' a>oo ' 








aoi 


= A 


aoi 




aio 




O10 




. an . 




. an . 



(27) 



The coefficients bj k introduced in Eq. (|27l) are the overlap 
integrals between the functions do(E) and d\(E): 



Ojk 



dEd*{E)d k {E). (28) 
hi = 1 due to the 



It holds that boi — b\ an d boo 
normalization. 

The system of algebraic equations (|27]l has a nontrivial 
solution provided that the eigenvalues A are solutions of 
the secular equation: 



A 



A 3 + VX-V 2 =0, 



where 



= 600O11 - ^oioio- 



(29) 



(30) 



The fourth-order polynomial in Eq. (1291) can be written as 
a product of the second-order polynomials (A 2 —T>){X 2 — 
X + V). This allows to find its roots: 



A x , 2 = ±VV, 



^3,4 



± 



V. 



(31) 



As the negativity TV is given by the amount of negative- 
ness in the eigenvalues A, we have 



N = VV. 



(32) 



Alternative and more intuitive derivation of the formula 
in Eq. (|32j) can be found in Appendix A invoking the 
decomposition of functions do(E) and di(E). 

In parallel to the entanglement, quantum discord [23[ 
has been discussed in the last years for systems composed 
of several parts [UJ • Discord quantifies the amount of in- 
formation in the whole system that cannot be extracted 
using quantum measurements at separated parts. Pro- 
vided that a bipartite system is in a pure state quantum 
discord is quantified by entropy S of entanglement. The 
entropy S of entanglement is given by the entropy of re- 
duced statistical operator g a of atom a that takes the 
form 



boo bio 
boi bn 



(33) 



exploiting the coefficients bj k . The eigenvalues A3,4 writ- 
ten in Eq. (131 j) naturally give also the eigenvalues of ma- 
trix g a and so they can be conveniently used in expressing 
the entropy S. The entropy S of entanglement is given 
by the usual formula S = — J2j=3 4 A? l°S2(^j)j 1°S2 be- 
ing the logarithm of base two. This formula provides us 
the following expression: 



-4 



log 2 (2?) + Vl-42?log 2 



1 + VI - 4X> 



1 - VI -42? 



(34) 



Here, determinant T> of the matrix g a is given in Eq. (|30l) . 

Combining Eqs. (|3"2"|) and (IM)) , the entropy S of entan- 
glement can be expressed as a monotonous function of 
negativity N (see Fig. [5]): 



S = -log, (AT) log, ; . 

2V ' 2 52 \1- V1-4AT 2 J 

(35) 

The curve in Fig. [2] reveals that both quantities can be 
equally well used for the quantification of entanglement 
in the considered system. 

The negativity can also be expressed in terms of 
eigenvalues of the Schmidt decomposition of the state 
in the long-time limit. Substituting Eq. (|3TJ)) into 
Eq. (|32p . we arrive at the useful formula for negativity 
N: 



N 



\ 



I 1 

9 Yl i b 3jhk - b jk b kj }. 



(36) 



j,k=0 



Further substitution for the coefficients bj k from Eq 
provides the negativity N depending on the reduced sta- 
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FIG. 2: Entropy S of entanglement as a function of negativity 
N in the interval [0, 0.5] of attainable values of TV. 



tistical operator g b of the continuum: 



N 



dE / dE'\g b {E,E 



\E,E>) = J2 dk(E)dt(E'). 



,(37) 
(38) 



fc=0,l 



Using the coefficients an d \AI of the Schmidt de- 
composition of the state \tp), the formula (|3"T|) can be 
recast into the simple form: 



N= ^hX 4 



(39) 



The formula for negativity N in Eq. (f3"2"j) can even be 
used for finite times t, in which discrete states of atom b 
are populated. In this case, the formula in Eq. (|28]l has 
to be replaced by the more general one: 



&i*=X)c;,c w + J dEd*(E)d k (E). 



(40) 



IV. QUADRATIC NEGATIVITY AND ITS 
SPECTRAL DENSITY 

The substitution of expression in Eq. (|3T))) into the for- 
mula (|32|) for negativity N gives us an expression that 
indicates the existence of quadratic negativity N q as a 
measure of entanglement that allows to introduce a spec- 
tral density [25| : 



AT, = 4iV 2 = 4(& 00 &ii - &oiM- 



(41) 



The use of expressions (|28|) for the coefficients bjk allows 
us to rewrite the formula in Eq. (|4l"T) as: 



n q = 2 J dEg{E) J dE' g(E')n q (E, E'), (42) 

where g(E) gives the density of states \E) in the contin- 
uum: 



The spectral density n q (E,E') of quadratic negativity 
introduced in Eq. (|42j) is obtained in the form: 



n q (E,E') = 



1 



q(E)q(E>) 



]T \d 3 {E)\ 2 \d k {E') 



j,k=0 



j,k=0 



d*(E)d k (E)d* k (E')d J (E') 



(44) 



The value of spectral density n q (E 7 E') of quadratic neg- 
ativity gives the value of quadratic negativity of a qubit- 
qubit system composed of the states {|0) a ,|l) Q } and 
{\E),\E'}}. According to Eq. (g2J), the quadratic neg- 
ativity N q is given as a weighted sum of quadratic qubit- 
qubit negativities between the two-level atom a and all 
possible qubits embedded inside the continuum. This in- 
terpretation is important from the physical point of view, 
because it allows to interpret the overall entanglement as 
composed of individual spectral contributions. We note 
that values of both the quadratic negativity N q and its 
density n q lie in the interval [0, 1]. We also note that an 
alternative normalization in the definition (|44[) of den- 
sity n q of quadratic negativity is possible. It is based 
on substituting the factor l/[g{E)g{E')\ by the factor 
4/[g(E) + g(E')] 2 . However, this 'mathematically more 
compact' normalization is not suitable for indicating en- 
tanglement in the case of qubits with considerably differ- 
ent values of the probability densities g(E) and g(E'). 

Experimental determination of the density n q of 
quadratic negativity has to take into account a finite 
resolution AE of frequencies of free electrons. That 
is why, it is convenient to introduce a series of exper- 



imental quadratic negativities N, 



0) 



1,2,..., that 



are obtained after spectral filtering of a free electron by 
using i filters positioned at the central frequencies Ek, 
k = 1, . . . ,i: 



N^(E 1 ,...,E i ) = 
b$(E 1 ,...,E i )b ( $(E 1 ,.. 



E i )-\b$(E 1 ,...,E i 



b$(E 1 ,...,E i )+b?>{E 1 ,...,E,, 



nl/2 



(45) 



The coefficients b^ {E\, . . . , Ei) occurring in Eq. (|45|) de- 
pend on the experimental frequency width AE and are 
given as: 

b®{E u ...,£,)=£ / dE'd*(E')d k (E'). 

l=1 JEi-AE/2 

(46) 

We note that the last term in the expression (|45j) origi- 
nates in the normalization of the considered state. 



ENTANGLEMENT GENERATION 



g(E) = J2\dAE)f 

3=0 



(43) The entanglement between electrons at atoms a and b 

is generated by the dipole-dipole interaction that is char- 
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acterized by the coefficients J a b and J. This means that 
two different channels of the entanglement generation ex- 
ist. In the first channel, the entanglement among the 
discrete states at atoms a and b is formed due to the 
dipole-dipole interaction described by the coefficient J a b 
first. Subsequently, this entanglement is transferred to 
the continuum of states \E) using either the Coulomb 
interaction (V) or the optical dipole interaction (/ita^). 
The second channel is based on the dipole-dipole interac- 
tion (J) between the excited discrete state |l) a at atom 
a and the continuum of states \E) at the ionization atom 
b. 

The dynamics of the system is such that an electron 
at atom b gradually 'leaks' into the continuum of states 
\E). The probability of finding this electron in a combi- 
nation of discrete states |0)& and 1 1) t, decreases roughly 
exponentially. After a sufficiently long time, this proba- 
bility is practically zero, the electron is fully ionized and 
its long-time spectrum completely characterizes its state. 
On the other hand, the electron at atom b periodically os- 
cillates between its discrete states in a stationary optical 
field at the Rabi frequency. The entanglement between 
the bound electron at atom a and the ionized electron at 
atom b is formed during the period of ionization and is 
'frozen' as soon as atom b is completely ionized. At this 
instant, the entanglement reaches its long-time limit, but 
superimposed periodic oscillations are possible under cer- 
tain conditions (see below). 

Let us concentrate on the first channel. Both electrons 
at atoms a and b being initially in their ground states 
gradually move into their excited states |l) a |0)b, |0) a |l)f>, 
and |l) a |l)(, due to the interaction with the stationary 
optical field [see Fig. |3Ja)]. The entanglement between 
discrete states arises from the dipole-dipole interaction 
between the states |l) |Q)b and |0) o |l)(,. The probabili- 
ties |cio| 2 and |coi| 2 affiliated to these states periodically 
return to zero with a period that decreases with the in- 
creasing values of \ J a b\, \Ha,otL\, and |/Z(,a:i,|. At these in- 
stants, highly entangled states occur and their quadratic 
negativities Nq quantifying entanglement among discrete 
states reach local maxima [see Fig. Efb)]. Provided that 
the probabilities of the ground state |0) a |0)(, and the state 
with both electrons excited are balanced ( | coo 1 2 ~ l c n| 2 )j 
the quadratic negativity Njf reaches its maximum value 
one. The quadratic negativity Nq oscillates between its 
maximum and zero during the time evolution. The en- 
tanglement between the discrete states at atom a and the 
continuum of states at atom b arises as a consequence 
of the interaction of the continuum of states \E) with 
the discrete states |0)& and The quadratic nega- 

tivity N q appropriate for this entanglement typically in- 
creases during the time evolution and gradually reaches 
its long-time value, as documented in Fig.[3Jb). However, 
weak oscillations may occur in this evolution. The overall 
quadratic negativity N$ that characterizes the entangle- 
ment between atoms a and b including all states, behaves 
similarly as the quadratic negativity N q comprising only 
the continuum of states. As a rule of thumb, a slightly 



(a) 





(b) 



FIG. 3: Temporal evolution of (a) probabilities | coo | 2 (solid 
curve with *), |cio| 2 (solid curve with o), |coi| 2 (solid curve 
with A), |cn| 2 (solid curve with o) of detecting two elec- 
trons in the appropriate discrete states and their sum E 
(E = X)j=o,i 2fe=o,i l c J fc | 2 ' solid curve ) and ( b ) quadratic 
negativities N q (solid curve with *), Nq (solid curve with A), 
and Nq (solid curve); /i n ai = \i b a.L = J a b = V = 0.05, 
H = J = 0,E a = E b =E L = l. 



stronger optical pumping of atom a compared to atom 
b (/j, a > /it) results in greater values of the long-time 
quadratic negativity Ng. 

In the second channel, the entanglement is generated 
directly by the dipole-dipole interaction between the ex- 
cited state |1)„ and the continuum of states \E). The 
ability to generate the entanglement is weaker compared 
to the first channel. 'Transfer of entanglement' can be ob- 
served also here and so nonzero values of the quadratic 
negativity N~ are found during the temporal evolution 
[see Fig. [4]. Even the maximum entangled discrete states 
(Nq = 1) can be reached. This clearly shows that there 
exists a strong 'back-action' from the 'reservoir' contin- 
uum of states \E) towards the discrete states \0) b and 
Otherwise, the observed temporal evolution is qual- 
itatively similar to that found in the first channel. 

Some general features of the behavior of quadratic neg- 
ativity N q in the long-time limit can be obtained even 
analytically. A detailed analysis of the long-time solu- 
tion in Eq. (p?Tj) has shown 11] that the coefficients &oo 
and &n giving the probabilities of finding an electron at 
atom a in the states |0) a and |l) a , respectively, can be 
expressed in the form: 



&oo(*) 

&i\(t) 



(1 



[6exp(i<5£t) + c.c] 
-a) — [bexp(iS^t) 



(47) 
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FIG. 4: Temporal evolution of quadratic negativities N q (solid 
curve with *), Ng (solid curve with A), and Nq (solid curve); 
J = 0.05, Jab = 0, values of the other parameters are the 
same as in the caption to Fig. [3] 



Constant a describes the steady-state parts of probabil- 
ities &oo and &n, whereas constant b gives the amount 
of probability that oscillates between the states |0) a and 
|1) Q at the Rabi frequency 5£. The symbol c.c. replaces a 
complex conjugate term. On the other hand, the cross- 
correlation coefficient i>oi can be written as 



boi{t) = c\ + c 2 exp(i6£t) + c 3 exp(-id^t); 



(48) 



ci, C2, and C3 being constants. 

Using the equality b 2 — — C2C3 valid in the model, 
we arrive at the following formula for the long-time 
quadratic negativity JVi*: 



JV*(t) = 4|a(l -«) - 2|6| 2 - 



cil 



C2 



M 



(1 — 2a)b — c\c2 



Cl b 2 

C2 



exp(z<5^t) 



(49) 



We can see from Eq. (|49[) that the quadratic negativity 
N^f is composed of a steady-state part and an oscillating 
part with the Rabi frequency 6£. However, the oscil- 
lating part is usually much smaller than the steady-state 
one. Even if atom a is resonantly pumped, the oscillating 
term in Eq. (1491) vanishes and we arrive at the simplified 
formula: 



N: 



It, res 



a(l-a)-2|6| 2 -| Cl | 2 -| C2 | 5 



\b\' 



C2 



(50) 

According to Eq. (f5*U)) . equal steady-state probabilities a 
and (1 — a) of detecting the electron at atom a in the 
states \0} a and |l) a , respectively, are needed to reach the 
maximum value of quadratic negativity Ng (a = 1/2). 
Moreover, nonzero values of constants |6|, |ci| and \c2\ 
lower the values of long-time quadratic negativity N q l . 

The numerical analysis of the long-time behavior of 
quadratic negativity Ng has revealed that the larger the 
values of dipole-dipole constants J a b and J are, the larger 
is the potential to generate highly entangled states. In 
order to arrive at high values of the quadratic negativity 
JVg*, the values of constants /mil and V have to be suffi- 
ciently small compared to the values of J and J a b- This 



can be physically explained as follows. The constants 
/ittL and V determine the speed of transfer of an electron 
at atom b into the continuum of states \E). If this speed 
is too fast, the electron at atom b has not enough time to 
create the entanglement with the electron at atom a. As 
a consequence, the entanglement between two electrons is 
weaker. This behavior is documented in Fig. [5] consider- 
ing both channels of entanglement generation. However, 
the graphs in Fig. [5] reveal that also greater values of the 
constants \ia l and V allow to reach strong entanglement 
under the condition /xa^ ~ V . The analysis of tempo- 
ral behavior of the system has shown that the movement 
of the electron at atom b into the continuum of states 
\E) is considerably slowed down in this case of balanced 
interactions /la^ and V. This slowing-down then gives 
enough time for the entanglement generation even for 
smaller values of the constants J a b and J. This regime is 
even preferred for the channel exploiting the constant J, 
as the graph in Fig. EJb) shows. 

Two channels based on the constants J a b and J mutu- 
ally 'interfere' in creating the entanglement between two 
electrons. This can be conveniently used for reaching 
greater values of the quadratic negativity Ng in regions, 
where the above described conditions are not met. Great 
values of the quadratic negativity Ng can be obtained in 
specific areas of the space spanned by the constants J a b 
and J, as illustrated in Fig. |6] 

We have considered the resonant pumping of atoms a 
and b up to now. The non-resonant pumping of both 
atoms makes the dynamics as well as the entanglement 
generation even more complex. Upon depending on con- 
ditions, the frequency detuning of atoms a and b may 
either support the entanglement creation or degrade it. 
A typical graph showing the behavior of quadratic nega- 
tivity N^ in dependence on the detunings A a and Ah is 
plotted in Fig. 



VI. SPECTRAL ENTANGLEMENT 

We illustrate typical properties of the spectral entan- 
glement considering the system characterized by param- 
eters mentioned in the caption to Fig. [3] In Fig. [51 the 
spectral density n q of quadratic negativity is plotted in 
the range of relative frequencies that covers two com- 
plex peaks occurring in the ionization spectrum (shown 
in Fig. in]). Strong spectral correlations inside the com- 
plex peaks as well as between different peaks are clearly 
visible. They mainly occur in spectral regions where the 
fast intensity variations occur (compare Figs. [8J and [9]). 

The experimental quadratic negativity iVq defined in 
Eq. (l4"5*j) represents the simplest experimentally accessible 
characteristics. As its definition indicates, the negativity 
./Vq 1 ^ depends on the experimental frequency resolution 
AE. It even holds that N^ 1] (E) -> for AE -> 0. This 
reflects the fact that at least a 'small' group of states \E) 
inside the frequency interval AE is needed to 'imprint' 
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FIG. 5: Topo graphs of long-time quadratic negativity Nq 
depending on optical pumping [ioil and strength V of the 
Coulomb interaction for (a) J ab = 0.001, J = 0, (b) J ab — 0, 
J = 0.001, and (c) J ab = J = 0.001; ^ a a L = /U&ai = 0.05, 
E a = Eb = El = 1. 



the entanglement. The wider the frequency interval AE 
is, the larger are the values of quadratic negativity Nq. 
As an example, the long-time 'distribution' of entangle- 
ment along the relative frequency axis (E — Eb)/T for 
the case studied in Fig. [3] is shown in Fig. [TDJ According 
to Fig. [10] there exist four spectral regions that consid- 
erably contribute to the formation of entanglement. If 
the frequency interval AE is sufficiently wide, the maxi- 
mum attainable values of quadratic negativity Nq' lt can 
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FIG. 6: Topo graphs of long-time quadratic negativity 

as it depends on dipole-dipole coupling constants J ab and J; 

fi a a L = fi b a L = i±a L = 0.05, V = 0.01, E a = E b — E L = 1. 
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FIG. 7: Topo graph of long-time quadratic negativity N l q x as 
a function of detunings A a and A b of atoms a and b, respec- 
tively; HaOL = /ifcCtL = 0.05, (ictL = 0.005, J ab = V = 0.001, 
J = 0, E L = 1. 





1.0 



FIG. 8: Topo graph of density of quadratic negativity 
showing qubit-qubit correlations in relative frequencies (E — 
E b )/Y and (E' — E b )/Y, values of parameters given in the 
caption to Fig. [3] are used. 
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FIG. 9: Long-time photoelectron ionization spectrum I • val- 
ues of parameters given in the caption to Fig. [3] are used. 
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(E-E„)/r 

FIG. 10: Long-time experimental quadratic negativity Nq 14 
as a function of relative frequency (E — Eb)/Y for AE/T = 
0.001 (solid curve), AE/T = 0.005 (solid curve with *), 
AE/T = 0.025 (solid curve with A), and AE /F = 0.05 (solid 
curve with o); r = 7r|V"| 2 + 7i~|J| 2 , values of parameters given 
in the caption to Fig. [3] are used. 



be approached. The comparison of the graph in Fig. [TU] 
with that in Fig. [9] giving the long-time photoelectron 
ionization spectrum J 1 * reveals that two spectral regions 
in the middle are crucial for constituting the entangle- 
ment between two electrons. 

We note that the experimental quadratic negativity 
Nq 1 * 1 is time-independent in the long-time limit provided 
that atom a is resonantly pumped. We remind that this 
is not the case of conditional long-time photoelectron ion- 
ization spectra Jg* and ij* obtained for atom a being in 
the ground (|0) a ) and the excited (|l) a ) state, respec- 



tively (for details, see [11 



The spectral correlations of entanglement as theoreti- 
cally described by the density n q (E, E') of quadratic neg- 
ativity can be experimentally revealed measuring the ex- 

perimental quadratic negativity Nq (E,E') introduced 
in Eq. (|45|) . As the considered example documents in 
Fig. QTJ two kinds of the spectral correlations of entan- 
glement may be distinguished. Strong correlations are 
found among the frequencies E and E' lying inside one 
spectral peak, but different sub-peaks [see Fig. [TOT a)]. 
On the other hand, strong correlations occur also for the 
frequencies E and E' localized inside the neighbor spec- 
tral peaks. Here, the correlations are observed inside the 
lower-frequency sub-peaks of two neighbor spectral peaks 
as well as inside the upper-frequency sub-peaks of the 
neighbor peaks [see Fig. fTUf b)]. This example illustrates 
richness of the internal spectral structure of entangled 
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(b) 

FIG. 11: Topo graphs of long-time experimental quadratic 
negativities iV, 2 '' 1 * depending on relative frequencies (E — 
Eb)/Y and (E' — Et)/T and showing the correlations (a) in- 
side one spectral region and (b) between two different spec- 
tral regions; AE/T = 0.001, values of the other parameters 
are written in the caption to Fig. [3] 



stated in the investigated system. 



VII. CONCLUSIONS 

The entanglement between two electrons in an auto- 
ionization atom and a neighbor two-level atom has been 
investigated. An expression for the negativity of a bi- 
partite system composed of a qubit and a general system 
including both the discrete and continuum levels has been 
derived. The spectral density of quadratic negativity has 
been introduced to study the spectral features of entan- 
glement. It has allowed to decompose the overall entan- 
glement into the qubit-qubit entanglement of the consti- 
tuting parts. Also the concept of experimental quadratic 
negativities has been introduced. It has been shown 
that the dipole-dipole interaction creates the entangle- 
ment between electrons until one of them is completely 
ionized. This puts restrictions to the strength of ion- 
ization paths in the auto-ionization atom. However, the 
balancing of two ionization paths in the auto-ionization 
atom results in a lower ionization speed that is in favor 
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of the entanglement generation. Highly entangled states 
stable for long times are then reached. The entangle- 
ment is spectrally 'concentrated' below the peaks of the 
long-time ionization spectra. Strong correlations have 
been found for pairs of frequencies localized inside one 
spectral peak as well as when two frequencies have been 
positioned below the neighbor peaks. 



Appendix A: Alternative derivation of the formula 
(f32|) for negativity TV 

We may conveniently decompose the functions d (E) 
and d\ (E) characterizing an ionized electron at atom b in 
a suitable orthonormal basis formed by functions fo(E) 
and fi(E). In this basis, the problem of quantifying en- 
tanglement between the two-level system a and the sys- 
tem b with the continuum of states is reduced to the 
problem of quantifying the entanglement in a qubit-qubit 
system. The appropriate basis functions fo(E) and fi{E) 
can be constructed along the following recipe: 



ME) 
fi(E) 



dy{E) 

Vb~oo~ ' 
-b w d (E) 



b 00 di(E) 



froobii ^ \b<n | 



(Al) 



the coefficients bjk have been defined in Eq. (1281) . The 
inverse transformation to that wr itten in Eq. dATj) can 
be derived in the form: 



d (E) 
di(E) 



&00 — V^OOj CflO 

bio] 



\bio\ 2 )/b 00 . 



aoofo{E), 

a w fo(E) + a n h{E); (A2) 
bio/Vkm, and an = {b 00 bu - 



Using new basis vectors \0))b and \l))b in the contin- 
uum of states at atom b, 



dEfj(E)\E), j — 0, 1, 



(A3) 



the state vector in Eq. 
following long-time form: 



can be recast into the 



\ip) U = a o|0) |0)) 6 + a 10 |l)„|0» 6 + a 1 i|l) a |l)) 6 . 

(A4) 

The state vector \ip) lt can be considered as a state of 
two qubits, a and b. The partially transposed statistical 
operator g Ta , transposed with respect to the indices of 
atom a, can be written in the following matrix form: 



Q Ta = 



oo 



aooaio aooaii 




laiol 2 



aioan 



(A5) 



The secular equation for the matrix g Ta can be obtained 
in the form (A 2 - V) (A 2 - p\ + V) = 0, where p 



and V has been defined in Eq. ([313} . 



The 
-y/V, 



only negative solution of the secular equation, A 
gives the formula for negativity N given in Eq. 

Acknowledgments 



Support by projects COST OC 09026, 1M06002 
and Operational Program Research and Development 
for Innovations - European Social Fund (project 
CZ. 1.05/2. 1.00/03.0058) of the Ministry of Education of 
the Czech Republic are acknowledged. Also support by 
project PrF-201 1-009 of Palacky University is acknowl- 
edged. 



[1] J. Matulewski, A. Raczyhski, and J. Zaremba, Phys. Rev. 

A 68, 013408 (2003). 
[2] G. S. Agarwal, S. L. Haan, and J. Cooper, Phys. Rev. A 

29, 2552 (1984). 
[3] W. Leohski, R. Tanas, and S. Kielich, J. Opt. Soc. Am. 

B 4, 72 (1987). 
[4] W. Leoriski and R. Tanas, J. Phys. B: At. Mol. Opt. 

Phys. 21, 2835 (1988). 
[5] W. Leoriski, R. Tanas, and S. Kielich, J. Phys. D: Appl. 

Phys. 21, S125 (1988). 
[6] U. Fano, Phys. Rev. 124, 1866 (1961). 
[7] K. Rzazewski and J. H. Eberly, Phys. Rev. Lett. 47, 408 

(1981). 

[8] P. Lambropoulos and P. Zoller, Phys. Rev. A 24, 379 
(1981). 

[9] A. Luks, V. Pefinova, J. Pefina Jr., J. Kfepelka, and 
W. Leoriski, in Wave and Quantum Aspects of Contem- 
porary Optics: Proceedings of SPIE, Vol. 7746, edited 
by J. Miillerova, D. Senderakova, and S. Jurecka (SPIE, 
Bellingham, 2010), p. 77460W. 



[10] J. Pefina Jr., A. Luks, W. Leoriski, and V. Pefinova, 

Phys. Rev. A 83, 053416 (2011). 
[11] J. Pefina Jr., A. Luks, W. Leoriski, and V. Pefinova, 

Phys. Rev. A 83, 053430 (2011). 
[12] J. Pefina Jr., A. Luks, V. Pefinova, and W. Leoriski, Opt. 

Express 18, 17133 (2011). 
[13] J. Pefina Jr., A. Luks, V. Pefinova, and W. Leoriski, J. 

Russian Laser Res. 33, 212 (2011). 
[14] L. Journel, B. Rouvellou, D. Cubaynes, J. M. Bizau, F. J. 

Willeumier, M. Richter, P. Sladeczek, K.-H. Selbman, 

P. Zimmerman, and H. Bergerow, J. de Physique IV 3, 

217 (1993). 

[15] A. Raczyriski, M. Rzepecka, J. Zaremba, and S. Zieliriska- 
Kaniasty, Optics Communications 266, 552 (2006). 

[16] W. Leoriski, J. Opt. Soc. Am. B 10, 244 (1993). 

[17] E. A. Silinsh and V. Capek, Organic Molecular Crys- 
tals: Interaction, Localization and Transport Phenomena 
(Oxford University Press/ American Institute of Physics, 
1994). 

[18] D. Bouwmeester, A. Ekert, and A. Zeilinger, The Physics 



11 



of Quantum Information (Springer, Berlin, 2000). 

[19] M. A. Nielsen and I. L. Chuang, Quantum Computa- 
tion and Quantum Information (Cambridge Univ. Press, 
Cambridge, 2000). 

[20] F. R. Gantmacher, The Theory of Matrices 
(AMS Chelsea publishing, Providence, 2000). 

[21] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78, 5022 
(1997). 

[22] G. Vidal and R. F. Werner, Phys. Rev. A 65, 032314 



(2002). 

[23] H. Ollivier and W. H. Zurek, Phys. Rev. Lett. 88, 017901 
(2001). 

[24] A. Al-Qasimi and D. F. V. James, Phys. Rev. A 83, 
032101 (2011). 

[25] W. Diir, J. I. Cirac, and R. Tarrach, Phys. Rev. Lett. 83, 
3562 (1999). 



